Unexpected high yield of acrolein underlies the importance of the hydrogen-abstraction mechanism in photooxidation of allyl methyl sulfide (AMS).

This work explores theoretically the gas phase oxidation of allyl methyl sulfide (AMS, H2CCHCH2SCH3) initiated by •OH radicals, focusing on the H-abstraction pathway at the M06-2X-D3/aug-cc-pVTZ and MN15/aug-cc-pVTZ levels of theory (m06Tz and mn15Tz). The formation of a prereactive complex (PRC) is involved in H-abstraction processes with two potential directions of approach for the OH radical, denoted as "α" and "ß". The PRCs, demonstrate increased reactivity, primarily due to the interaction between the sulfur atoms and the hydroxyl hydrogen. A scheme for the H-abstraction mechanism that supports the experimentally identified products and predicts the formation of some S-containing low volatility products is proposed. The comparison of the potential energy surface (PES) between the double bond addition and H-abstraction paths in the AMS molecule shows that at the m06Tz level of theory, the H-abstraction on C3 and the addition to C1 have nearly the same profile of energy, while at the mn15Tz level, the minimum energy channel is the addition to C1. The theoretical rate coefficient for each reaction channel was calculated, considering the formation of a PRC prior to reaching the transition state of each channel and assuming thermal equilibrium between reactants and the PRC. The rate constants were calculated in a multi-TS/multi-conformer way at the SVECV-f12/m06Tz and SVECV-f12/mn15Tz levels of theory. The SVECV-f12 method is consistent in its predictions in both systems and exhibits only minor deviations from the experimental rate constants. Despite some specific differences due to the DFT method supporting the SVECV-f12 calculations, both methodologies predict a significant H-abstraction contribution in the AMS + OH gas phase reaction, which explains the high formation yield for acrolein determined experimentally.

A new chapter in the never ending story of cycloadditions: The puzzling case of SO2 and acetylene.

A comprehensive study of the different classes of cycloaddition reactions ([3+2], [2+2], and [2+1]) of SO2 to acetylene and ethylene has been performed using density functional theory (DFT) and composite wavefunction methods. The [3+2] cycloaddition reaction, that was previously explored in the context of the cycloaddition of thioformaldehyde S-methylide (TSM) to ethylene and acetylene, proceeds in a concerted way to the formation of stable heterocycles. In this paper, we extend our study to the [2+2] and [2+1] cycloadditions of SO2 to acetylene, which would produce 1,1-oxathiete-2-oxide and thiirene-1,1-dioxide, respectively. One of the main conclusions is that cyclic 1,1-oxathiete-2-oxide can open through a relatively easy breaking of the SO single bond and rearrange toward sulfinyl acetaldehyde (SA). The SA molecule can easily undergo several internal rearrangements, which eventually lead to sulfenic acid and sulfoxide derivatives of ethenone, 1,2,3-dioxathiole, and CO plus sulfinylmethane. The most probable path, however, produces 2-thioxoacetic acid, whose derivatives (or those of the corresponding acetate) are usually obtained by Willgerodt-Kindler-type sulfuration of acetates. This product can in turn decompose, leading to the final products CO2 and H2CS. Comparison of this decomposition path with that of 2-amino-2-thioxoacetic acid shows that the process occurs through different H-transfer processes.

Atmospheric oxidation of furanones by •OH and •Cl radicals: In situ FTIR rate coefficient determinations, SAR and theoretical studies.

Gas-phase kinetics of the overall reactions of •OH and •Cl radicals with dihydrofuran-3(2H)-one (oxolan-3-one) and dihydro-2-methyl-3(2H)-furanone (2MTHF-3-one) were studied at 298 K and atmospheric pressure. The rate coefficients were determined using the relative method in a 480 L multipass glass reactor coupled to an FT-IR detection system. The rate coefficients found for oxolan 3-one and 2MTHF-one with •OH radicals (k1 and k2) and with •Cl atoms (k3 and k4) at 298 K and atmospheric pressure (in cm3 molecule-1 s-1) were: k1 = (1.86 ± 0.29) × 10-11, k2 = (2.64 ± 0.47) × 10-11, k3= (1.15 ± 0.28) × 10 -10, and k4 = (1.33 ± 0.32) × 10-10, respectively. Reactivity trends were developed by comparison with other similar structures and Fukui indices employed to determine the reactivity of different sites on the ring. The singularity of the reaction with •OH was assessed by computational studies which showed the formation of several stable hydrogen bonded complexes, explaining the difference with the reaction with the •Cl atom. SAR estimations of the rate coefficients were calculated and compared to the experimental values.

Experimental and Theoretical Kinetic Studies of the Ozonolysis of Selected Allyl Sulfides (H2C═CHCH2SR, R = CH3, CH3CH2): The Effect of Nascent OH Radicals.

Rate coefficients of the O3-initiated oxidation of allyl methyl sulfide (H2C═CHCH2SCH3, AMS) and allyl ethyl sulfide (H2C═CHCH2SCH2CH3, AES) were determined at atmospheric conditions by "in situ" FTIR. The relative kinetic experiments were performed using methylcyclohexane (McH) and carbon monoxide (CO) as nascent OH radical scavengers and in the absence of any scavenger, to determine the impact that the formation of OH radicals has on the rate coefficients. In the absence of scavengers, values of kAMS+O3 = (5.23 ± 3.57) × 10-18 and kAES+O3 = (5.76 ± 1.80) × 10-18 cm3 molecule-1 s-1 were obtained. In the presence of the scavengers, however, the rates decreased to kAMS+O3+McH = (3.92 ± 1.92) × 10-18 and kAMS+O3+CO = (2.63 ± 0.47) × 10-18 cm3 molecule-1 s-1 for AMS, and kAES+O3+McH = (4.78 ± 1.38) × 10-18 and kAES+O3+CO = (3.50 ± 0.27) × 10-18 cm3 molecule-1 s-1 for AES. On the basis of these results, we have decided to recommend the values obtained using CO as scavenger as those best representing the rate coefficient for the reactions of O3 with AMS and AES. The reaction mechanism was explored using DFT and post-Hartree-Fock computational methods. It is shown that the barrier for the common cyclization to primary ozonide (-3.7 ± 0.1 kcal mol-1) followed by other reactions, as well as that for the reaction with the sulfur atom (-5.1 ± 0.1 kcal mol-1), is small and quite close, meaning that both reaction paths should contribute significantly to the global reaction rate.

Correcting the Experimental Enthalpies of Formation of Some Members of the Biologically Significant Sulfenic Acids Family.

Sulfenic acids are important intermediates in the oxidation of cysteine thiol groups in proteins by reactive oxygen species. The mechanism is influenced heavily by the presence of polar groups, other thiol groups, and solvent, all of which determines the need to compute precisely the energies involved in the process. Surprisingly, very scarce experimental information exists about a very basic property of sulfenic acids, the enthalpies of formation. In this Article, we use high level quantum chemical methods to derive the enthalpy of formation at 298.15 K of methane-, ethene-, ethyne-, and benzenesulfenic acids, the only ones for which some experimental information exists. The methods employed were tested against well-known experimental data of related species and extensive CCSD(T) calculations. Our best results consistently point out to a much lower enthalpy of formation of methanesulfenic acid, CH3SOH (ΔfH0(298.15K) = -35.1 ± 0.4 kcal mol-1), than the one reported in the NIST thermochemical data tables. The enthalpies of formation derived for ethynesulfenic acid, HC≡CSOH, +32.9 ± 1.0 kcal/mol, and benzenesulfenic acid, C6H5SOH, -2.6 ± 0.6 kcal mol-1, also differ markedly from the experimental values, while the enthalpy of formation of ethenesulfenic acid CH2CHSOH, not available experimentally, was calculated as -11.2 ± 0.7 kcal mol-1.

Dipolar 1,3-cycloaddition of thioformaldehyde S-methylide (CH2 SCH2 ) to ethylene and acetylene. A comparison with (valence) isoelectronic O3 , SO2 , CH2 OO and CH2 SO.

Methods rooted in the density functional theory and in the coupled cluster ansatz were employed to investigate the cycloaddition reactions to ethylene and acetylene of 1,3-dipolar species including ozone and the derivatives issued from replacement of the central oxygen atom by the valence-isoelectronic sulfur atom, and/or of one or both terminal oxygen atoms by the isoelectronic CH2 group. This gives rise to five different 1,3-dipolar compounds, namely ozone itself (O3 ), sulfur dioxide (SO2 ), the simplest Criegee intermediate (CH2 OO), sulfine (CH2 SO), and thioformaldehyde S-methylide (CH2 SCH2 , TSM). The experimental and accurate theoretical data available for some of those molecules were employed to assess the accuracy of two last-generation composite methods employing conventional or explicitly correlated post-Hartree-Fock contributions (jun-Cheap and SVECV-f12, respectively), which were then applied to investigate the reactivity of TSM. The energy barriers provided by both composite methods are very close (the average values for the two composite methods are 7.1 and 8.3 kcal mol-1 for the addition to ethylene and acetylene, respectively) and comparable to those ruling the corresponding additions of ozone (4.0 and 7.7 kcal mol-1 , respectively). These and other evidences strongly suggest that, at least in the case of cycloadditions, the reactivity of TSM is similar to that of O3 and very different from that of SO2 .

T- and pH-dependent OH radical reaction kinetics with glycine, alanine, serine, and threonine in the aqueous phase.

Glycine, alanine, serine, and threonine are essential amino acids originating from biological activities. These substances can be emitted into the atmosphere directly. In the present study, the aqueous phase reaction kinetics of hydroxyl radicals (˙OH) with the four amino acids is investigated using the competition kinetics method under controlled temperature and pH conditions. The following T-dependent Arrhenius expressions are derived for the ˙OH reactions with glycine, k(T, H2A+) = (9.1 ± 0.3) × 109 × exp[(-2360 ± 230 K)/T], k(T, HA±) = (1.3 ± 0.1) × 1010 × exp[(-2040 ± 240 K)/T]; alanine, k(T, H2A+) = (1.4 ± 0.1) × 109 × exp[(-1120 ± 320 K)/T], k(T, HA±) = (5.5 ± 0.2) × 109 × exp[(-1300 ± 200 K)/T]; serine, k(T, H2A+) = (1.1 ± 0.1) × 109 × exp[(-470 ± 150 K)/T], k(T, HA±) = (3.9 ± 0.1) × 109 × exp[(-720 ± 130 K)/T]; and threonine, k(T, H2A+) = (5.0 ± 0.1) × 1010 × exp[(-1500 ± 100 K)/T], k(T, HA±) = (3.3 ± 0.1) × 1010 × exp[(-1320 ± 90 K)/T] (in units of L mol-1 s-1). The energy barriers of the ˙OH-induced H atom abstractions were simulated by the density functional theory (DFT) calculation performed with GAUSSIAN using the method of M06-2X and the basis set of 6-311++G(3df,2p). According to the calculation results, the -COOH and -NH3+ groups with strong negative inductive effects increase the energy barriers and thus decrease the ˙OH reaction rate constants. In contrast, the presence of a -OH or -CH3 group with weak negative or positive inductive effects can reduce energy barriers and hence increase the ˙OH reaction rate constants. To improve the previous structure-activity relationship, the contribution factors of -NH3+ at Cα-atom and Cß-atom are determined as 0.07 and 0.15, respectively. Aqueous phase ˙OH oxidation acts as an important sink of the amino acids in the atmosphere, and can be accurately described by the obtained Arrhenius expressions under atmospheric conditions.

Degradation mechanism of 2-fluoropropene by Cl atoms: experimental and theoretical products distribution studies.

The gas-phase reaction products of 2-fluoropropene (2FP) with Cl atoms have been determined for the first time at 298 K and atmospheric pressure using a 1080 L quartz-glass photoreactor coupled with in situ FTIR spectroscopy to monitor reactants and products. Acetyl fluoride and formyl chloride were observed as the main products with yields of (106 ± 10)% and (100 ± 11)%, respectively. Electronic structure calculations of reactants, intermediates, products and transition states on a detailed mechanism of the reaction were performed by DFT procedures (BMK, M06, M062X/D3), as well as accurate composite methods on both the addition and abstraction reaction channels. From the joint experimental and theoretical studies, we concluded that the reaction occurs primarily via addition to the Cα carbon, with a smaller participation of the addition on the Cß carbon, which is not produced directly from the separated reactants but from the CH3CFCH2Cl intermediate radical through a submerged transition state. The abstraction channel occurs at larger energies than the addition ones, and also presents a submerged transition state, with a lower barrier. No products arising from this channel are expected. The proposed mechanism explains also why formaldehyde, predicted as a product by former theoretical studies, is not found among the experimental products. The atmospheric implications of the reaction products are assessed.

Isomerization and Fragmentation Reactions on the [C2SH4] Potential Energy Surface: The Metastable Thione S-Methylide Isomer.

Thione S-methylide, parent species of the thiocarbonyl ylide family, is a 1,3-dipolar species on the [C2SH4] potential energy surface, not so much studied as its isomers, thiirane, vinyl thiol, and thioacetaldehyde. The conrotatory ring-closure reaction toward thiirane was studied in the 90s, but no complete analysis of the potential energy surface is available. In this paper, we report a computational study of the reaction scheme linking all species. We employed several computational methods (density functional theory, CCSD(T) composite schemes, and CASSCF/CASPT2 multireference procedures) to find the best description of thione S-methylide, its isomers, and transition states. The barrier from thiirane to thione S-methylide amounts to 52.2 kcal mol-1 (against 17.6 kcal mol-1 for the direct one), explaining why thiocarbonyl ylides cannot be prepared from thiiranes. Conversion of thiirane to vinyl thiol implies a large barrier, supporting why the reaction has been observed only at high temperatures. Fragmentations of thiirane to S(3P) or S(1D) and ethylene as well as decomposition to hydrogen sulfide plus acetylene were also explored. Triplet and singlet open-shell species were identified as intermediates in the fragmentations, with energies lower than the transition state between thiirane and vinyl thiol, explaining the preference of the latter at low temperatures.

Reinvestigation of the Deceptively Simple Reaction of Toluene with OH and the Fate of the Benzyl Radical: The "Hidden" Routes to Cresols and Benzaldehyde.

In a previous work, we have investigated the initial steps of the reaction of toluene with the hydroxyl radical using several quantum chemical approaches including density functional and composite post-Hartree-Fock models. Comparison of H-abstraction from the methyl group and additions at different positions of the phenyl ring showed that the former reaction channel is favored at room temperature. This conclusion appears at first sight incompatible with the experimental observation of a lower abundance of the product obtained from abstraction (benzaldehyde) with respect to those originating from addition (cresols). Further reactions of the intermediate radicals with oxygen, water, and additional OH radicals are explored in this paper through theoretical calculations on more than 120 species on the corresponding potential energy surface. The study of the addition reactions, to obtain the cresols through hydroxy methylcyclodienyl intermediate radicals, showed that only in the case of o-cresol the reaction proceeds by addition of O2 to the ring, internal H-transfer, and hydroperoxyl abstraction and not through direct H-abstraction. For both p- and m-cresol, instead, the reaction occurs through a higher-energy direct H-abstraction, thus explaining in part the observed larger concentration of the ortho isomer in the final products. It was also found that the benzyl radical, formed by H-abstraction from the methyl group, is able to react further if additional OH is present. Two reaction paths leading to o-cresol, two leading to p-cresol, and one leading to m-cresol were determined. Moreover, in this situation, the benzyl radical is predicted to produce benzyl alcohol, as was found in some experiments. The commonly accepted route to benzaldehyde was found to be not the energetically favored one. Instead, a route leading to the benzoyl radical (and ultimately to benzoic acid) with the participation of one water molecule was clearly more favorable, both thermodynamically and kinetically.

Calculation of the Geometries and Infrared Spectra of the Stacked Cofactor Flavin Adenine Dinucleotide (FAD) as the Prerequisite for Studies of Light-Triggered Proton and Electron Transfer.

Flavin cofactors, like flavin adenine dinucleotide (FAD), are important electron shuttles in living systems. They catalyze a wide range of one- or two-electron redox reactions. Experimental investigations include UV-vis as well as infrared spectroscopy. FAD in aqueous solution exhibits a significantly shorter excited state lifetime than its analog, the flavin mononucleotide. This finding is explained by the presence of a "stacked" FAD conformation, in which isoalloxazine and adenine moieties form a π-complex. Stacking of the isoalloxazine and adenine rings should have an influence on the frequency of the vibrational modes. Density functional theory (DFT) studies of the closed form of FAD in microsolvation (explicit water) were used to reproduce the experimental infrared spectra, substantiating the prevalence of the stacked geometry of FAD in aqueous surroundings. It could be shown that the existence of the closed structure in FAD can be narrowed down to the presence of only a single water molecule between the third hydroxyl group (of the ribityl chain) and the N7 in the adenine ring of FAD.

Computational Evidence Suggests That 1-Chloroethanol May Be an Intermediate in the Thermal Decomposition of 2-Chloroethanol into Acetaldehyde and HCl.

The dehalogenation of 2-chloroethanol (2ClEtOH) in the gas phase with and without the participation of catalytic water molecules has been investigated using methods rooted into the density functional theory. The well-known HCl elimination leading to vinyl alcohol (VA) was compared to the alternative elimination route toward oxirane and shown to be kinetically and thermodynamically more favorable. However, the isomerization of VA to acetaldehyde in the gas phase, in the absence of water, was shown to be kinetically and thermodynamically less favorable than the recombination of VA and HCl to form the isomeric 1-chloroethanol (1ClEtOH) species. At the ωB97X-D/cc-pVTZ level of calculation, this species is more stable than 2ClEtOH by about 6 kcal mol-1 at 298 K, and the reaction barrier for VA to 1ClEtOH is 23 kcal mol-1 versus 55 kcal mol-1 for the direct transformation of VA to acetaldehyde. In a successive step, 1ClEtOH can decompose directly to acetaldehyde and HCl with a lower barrier (29 kcal mol-1) than that of VA to the same products (55 kcal mol-1). The calculations were repeated using a single ancillary water molecule (W) in the complexes 2ClEtOH_W and 1ClEtOH_W. The latter adduct is now more stable than 2ClEtOH_W by about 8 kcal mol-1 at 298 K, implying that the water molecule increased the already higher stability of 1ClEtOH in the gas phase. However, this catalytic water molecule lowers dramatically the barrier for the interconversion of VA to acetaldehyde (from 55 to 7 kcal mol-1). This barrier is now smaller than the one for the conversion to 1ClEtOH (which also decreases, but not so much, from 23 to 13 kcal mol-1). Thus, it is concluded that while 1ClEtOH may be a plausible intermediate in the gas phase dehalogenation of 2ClEtOH, it is unlikely that it plays a major role in water complexes (or, by inference, aqueous solution). It is also shown that neither in the gas phase nor in the cluster with one water molecule, the oxirane path is more favorable than the VA alcohol path. Additionally, a direct conversion of 2ClEtOH to 1ClEtOH through a transition state which resembles a VA molecule in a complex with a chlorine atom and a hydrogen atom on both sides of this planar species was found. This reaction path has also lower activation energy than the conversion to oxirane but not as low as the conversion to VA.

Theoretical study of the reactions of the hydroselenyl radical (HSe●) with the selenenic radical (HSeO●).

The formation of selenium species in some biological processes involves the generation of ionic and radical intermediates such as RSe●, RSe-, RSeO●, and RSeO-, among others. We performed a theoretical study of the possible mechanisms for the reaction of the two simplest Se radicals-the hydroselenyl (HSe●) and selenenic (HSeO●) radicals, in which the possible products, intermediates, and transition-state structures were investigated. Density functional theory (DFT) was applied at the B3LYP/6-311++G(3df,3pd) level and the Ahlrichs Coulomb fitting basis sets were employed with an effective core potential (ECP) for both Se atoms. The same procedure was used to calculate the electronic density. All calculations were also performed using the M06-2X functional, which describes weaker bonds better than B3LYP does. In the reaction of interest, the so-called CR complex (HSe····SeOH) is formed initially. After passing through the transition state TS1, cis-HSeSeOH is obtained as a product. If a low barrier is then overcome (passing through the transition state TS32), the trans-HSeSeOH species is obtained. The CR complex can also rearrange into the intermediate INT after overcoming the barrier presented by the transition state TS2. Additionally, the decomposition of INT to H2O and 1Se2 is possible through another transition state. This reaction is not included in this study. We also observed a second possible route for the conversion of INT to one of the HSeSeOH species; this route occurs through two pathways (with transition states TS31 and TS32). A comparison of some of the results with those obtained for sulfur analogs along the same pathways is also presented in this work. Graphical abstract Electronic envelopes for HSeO● and HSe● radicals.

Structural insights into human microsomal epoxide hydrolase by combined homology modeling, molecular dynamics simulations, and molecular docking calculations.

A new homology model of human microsomal epoxide hydrolase was derived based on multiple templates. The model obtained was fully evaluated, including MD simulations and ensemble-based docking, showing that the quality of the structure is better than that of only previously known model. Particularly, a catalytic triad was clearly identified, in agreement with the experimental information available. Analysis of intermediates in the enzymatic mechanism led to the identification of key residues for substrate binding, stereoselectivity, and intermediate stabilization during the reaction. In particular, we have confirmed the role of the oxyanion hole and the conserved motif (HGXP) in epoxide hydrolases, in excellent agreement with known experimental and computational data on similar systems. The model obtained is the first one that fully agrees with all the experimental observations on the system. Proteins 2017; 85:720-730. © 2016 Wiley Periodicals, Inc.

Theoretical insight into the mechanism for the inhibition of the cysteine protease cathepsin B by 1,2,4-thiadiazole derivatives.

Several cellular disorders have been related to the overexpression of the cysteine protease cathepsin B (CatB), such as rheumatic arthritis, muscular dystrophy, osteoporosis, Alzheimer's disease, and tumor metastasis. Therefore, inhibiting CatB may be a way to control unregulated cellular functions and prevent tissue malformations. The inhibitory action of 1,2,4-thiadiazole (TDZ) derivatives has been associated in the literature with their ability to form disulfide bridges with the catalytic cysteine of CatB. In this work, we present molecular modeling and docking studies of a series of eight 1,2,4-thiadiazole compounds. Substitutions at two positions (3 and 5) on the 1,2,4-thiadiazole ring were analyzed, and the docking scores were correlated to experimental data. A correlation was found with the sequence of scores of four related compounds with different substituents at position 5. No correlation was observed for changes at position 3. In addition, quantum chemistry calculations were performed on smaller molecular models to study the mechanism of inhibition of TDZ at the active site of CatB. All possible protonation states of the ligand and the active site residues were assessed. The tautomeric form in which the proton is located on N2 was identified as the species that has the structural and energetic characteristics that would allow the ring opening of 1,2,4-thiadiazole.

Improved homology model of cyclohexanone monooxygenase from Acinetobacter calcoaceticus based on multiple templates.

A new homology model of cyclohexanone monooxygenase (CHMO) from Acinetobacter calcoaceticus is derived based on multiple templates, and in particular the crystal structure of CHMO from Rhodococcus sp. The derived model was fully evaluated, showing that the quality of the new structure was improved over previous models. Critically, the nicotinamide cofactor is included in the model for the first time. Analysis of several molecular dynamics snapshots of intermediates in the enzymatic mechanism led to a description of key residues for cofactor binding and intermediate stabilization during the reaction, in particular Arg327 and the well known conserved motif (FxGxxxHxxxW) in Baeyer-Villiger monooxygenases, in excellent agreement with known experimental and computational data.

Mechanism of the organocatalyzed decarboxylative Knoevenagel-Doebner reaction. A theoretical study.

We have investigated important intermediates and key transition states of the organocatalyzed Knoevenagel condensation using density functional theory and two different basis sets (6-31 G(d,p) and 6-311++G(2df,2pd)), both in gas phase and simulating the bulk solvent (pyridine) using the PCM method. Calculated structures for reactants, intermediates, and key transition states suggest that the secondary amine catalyst is essential, both for activating the aldehyde for nucleophilic attack, and in the possible decarboxylation pathways. The calculated results are shown to agree with available experimental information. On the basis of the results obtained, the studied mechanism may be important in the understanding of vinylphenol production during malting and brewing of wheat and barley grains.

Quantum model of catalysis based on a mobile proton revealed by subatomic x-ray and neutron diffraction studies of h-aldose reductase.

We present results of combined studies of the enzyme human aldose reductase (h-AR, 36 kDa) using single-crystal x-ray data (0.66 A, 100K; 0.80 A, 15K; 1.75 A, 293K), neutron Laue data (2.2 A, 293K), and quantum mechanical modeling. These complementary techniques unveil the internal organization and mobility of the hydrogen bond network that defines the properties of the catalytic engine, explaining how this promiscuous enzyme overcomes the simultaneous requirements of efficiency and promiscuity offering a general mechanistic view for this class of enzymes.

In-silico nanobio-design. A new frontier in computational biology.

Nanobiology is a fast-emerging discipline that brings the tools of nanotechnology to the biological sciences. The introduction of new techniques may accelerate the development of highly specific biomedical treatments, increase their efficiency, and minimize their side effects. Introducing foreign bodies into the complex machinery of the human body is, however, a great and humbling challenge, as past experience has shown. In order for nanobiology to reach its full potential, we must devise a means to alter the properties of nanoparticles, as expressed in the human body, in a predictable manner. Computer-aided methods are the natural option to speed up the development of these technologies. Yet, the procedures for annotation and simulation of nanoparticle properties must be developed and their limitations understood before computational methods can be fully exploited. In this review we will compare the state of development of nanoscale simulations in the biological sciences to that of the computer-aided drug design efforts in the past, tracing a historical parallel between both disciplines. From this comparison, lessons can be learned and bottlenecks identified, helping to speed up the development of computer-aided nanobiodevice design tools.